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A solution of a hoary conundrum: the origin and properties of cosmic rays 

A. De Rujula a 

a Theory Division, CERN, 1211 Geneva 23, Switzerland 
Physics Department, Boston University, USA 

I discuss a theory of non-solar cosmic rays (CRs) based on a single type of CR source at all energies. All 
observed properties of CRs are predicted in terms of very simple and completely 'standard' physics. The source 
of CRs is extremely 'economical': it has only one parameter to be fitted to the enormous ensemble of all of the 
data. All other inputs are 'priors', that is theoretical or observational items of information independent of the 
properties of the source of CRs and chosen to lie in their pre-established ranges. 



1. Introduction 

A couple of years ago, I presented at a large 
cosmic-ray (CR) conference a preliminary version 
of the theory I am to summarize here. The cof- 
fee break subsequent to my talk was attended by 
hundreds of participants. A prominent CR theo- 
rist, who had heard my talk, told me: 'Your the- 
ory is unacceptable. It solves the whole problem'. 
He paused to make a large hand gesture embrac- 
ing the surrounding crowd, and added: 'What are 
then all these people going to do? Was he just jok- 
ing? Not that I am convinced that we have the 
complete solution to the problem. But, I shall ar- 
gue, we have a solution that Ockham would have 
favoured, for the reasons stated in the abstract. 

The standard theory posits that, up to the knee, 
CRs are produced by supernovae (SNe) [2], via 
shocks in the interaction of their roughly spher- 
ical non-relativistic ejecta with the interstellar 
medium (ISM) 0. The NASA website lists al- 
most 70000 refereed papers on 'cosmic ray(s)', 
most of them theoretical and 'standard'. As im- 
plied by a fraction of these papers, the standard 
theory has problems: it does not accelerate CRs 
up to the knee (e.g. 0]), the fraction of SN rem- 
nants perhaps compatible with the acceleration 
of nuclear CRs is insufficient to explain the CR 
luminosity of the Galaxy (e.g. 0), CRs produced 
by SN remnants (most of which are within the 
'solar circle') would diffuse outwards, generating 
a directional asymmetry which is not observed 
(e.g. [H]). CRs above the knee are mysterious. 



Our theory |2] is non-standard in that CRs are 
accelerated at all energies by the relativistic jets 
ejected by SNe. It is part of a 'unified theory of 
high-energy astrophysics' ^Q, based on the 'can- 
nonball' model [S] of the jets of accreting black 
holes and neutron stars, which also explains sim- 
ply the properties of gamma-ray bursts (GRBs) 
HJ, X-ray flashes (XRFs) their respective af- 
terglows (AGs) |11I12| . the gamma 'background' 
radiation |13|. the CR luminosity of our Galaxy 
|7I14| , and the properties of galaxy clusters har- 
bouring 'cooling flows' Many more and less 
self-referring citations are given in |7I9| . 

2. The Cannonball Model 

The 'cannon' of the CB model is analogous to 
the ones responsible for the ejecta of quasars and 
microquasars. Long-duration GRBs, for instance, 
are produced in ordinary core-collapse SNe by 
jets of CBs, made of ordinary-matter plasma, and 
travelling with high LFs, 7 ~ O(10 3 ). As a con- 
sequence of the initial star's rotation, an accre- 
tion disk is produced around the newly born com- 
pact object, either by stellar material originally 
close to the surface of the imploding core, or by 
more distant stellar matter falling back after the 
shock's passage |8I16) . A CB is emitted, as in mi- 
croquasars [TJ|, when part of the accretion disk 
falls abruptly onto the compact object. A sum- 
mary of the CB model is given in Fig. ^ 

Do SNe emit cannonballs? Until 2003 fig) , 
there was only one case with data good enough to 
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Figure 1. An 'artist's view' (not to scale) of 
the CB model of long-duration GRBs, XRFs and 
CRs. A core-collapse SN results in a compact ob- 
ject and a fast-rotating torus of non-ejected ma- 
terial. Matter (not shown) episodically accreting 
into the central object produces two narrowly col- 
limated beams of CBs; only some of the 'northern' 
CBs are depicted. As the CBs move through the 
'ambient light' surrounding the star, they forward 
Compton up-scatter its photons to GRB or XRF 
energies, depending on how close the line of sight 
is to the CBs' direction. Each CB produces a 
GRB 'pulse'. Later, a CB scatters ISM particles, 
which are isotropized by its inner magnetic field. 
In the SN rest system the particles arc boosted 
by the CB's motion: they have become CRs. The 
synchrotron radiation of the electrons is the late 
AG of the GRB or XRF. The CBs generate CRs 
all along their trajectories, in the galaxy, its halo, 
and beyond, as the CBs' collisions with the ISM 
slow them down. CRs are also forward-produced, 
diffusing thereafter in the local magnetic fields. 



tell: SN1987A, the core-collapse SN in the LMC, 
whose neutrino emission was seen. Speckle inter- 
ferometry data taken 30 and 38 days after the ex- 
plosion ^15] did show two back-to-back relativistic 
CBs. The approaching one was 'superluminal'. 

Are GRBs made by SNe? For long-duration 
GRBs, the answer is affirmative The first 

evidence for a GRB-SN association concerned 
SN1998bw, at redshift z = 0.0085, observed within 
the angular error towards GRB 980425. 

GRBs have 'afterglows': they are observable 
at radio to X-ray frequencies, for months after 
their 7-rays are seen. The optical luminosity of 
a 1998bw-like SN peaks at -15 (1+z) days. The 
SN light competes at that time with the AG of 
its GRB: it is not always detectable. It makes 
sense to test whether long-duration GRBs are 
associated with a 'standard torch' SN, akin to 
SN1998bw, 'transported' to their respective red- 
shifts. The test works optimally: for all — 15 
cases in which such a SN could be seen, it was 
seen (with varying significance) and for all — 15 
cases in which the SN could not be seen, it was 
not seen |llj . One could hardly do better. In 
practice SNe could not be observed at z > 1. 1. 

Naturally, truly 'standard torches' do not exist, 
but SN1998bw made such a good job of it that 
we could predict |20I21| the SN contribution to the 
AG in all six recent cases of early detection of the 
AGs of near-by GRBs. Besides 980425-1998bw, 
the most convincing association was provided by 
the spectroscopic discovery of a SN in the AG of 
GRBs 030329 for which even the night when the 
SN would be discovered was foretold [2] • 

In a CB-model analysis of GRBs and their AGs 
|9I11I12) we find that, within the pervasive cosmo- 
logical factor of a few, the long-GRB-SN associ- 
ation would be —1:1. Yet, current data are insuf- 
ficient to determine whether long-duration GRBs 
are associated with all core-collapse SNe (— 70% 
of all SNe, including Type II) or only with Type 
Ib/c SNe (-15% of core-collapse SNe). 

CB-model priors. The study of GRB AGs [TT1 
12 allowed us to extract, case by case, the initial 
LF, 70, and (less precisely) the initial mass and 
baryon number, M — iV B m p c 2 , of CBs, as well 
as 'environmental' quantities, such as the (highly 
varying) ISM density and the angle between the 
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jet of CBs and the observer. Typical values are: 

70 = E/(M c 2 ) « (10 3 ); iV B ~ O(10 50 ). (1) 

Mq is half the mass of Mercury, i.e. tiny, com- 
pared to the mass of the parent star. The 70 val- 
ues are roughly log-normally distributed around 
7o~1100 with a width at half-maximum extend- 
ing from 7^ ~ 630 to 7„ ~ 1400 [T2|. With these 
inputs, we could predict all properties of the indi- 
vidual 7 pulses of a GRB, each corresponding to 
a single CB whose electrons Compton-up-scatter 
the ambient light 0. The same inputs are used 
to predict the properties of CRs. 

While a CB exits from its parent SN and emits 
a GRB pulse, it is assumed [S] to be expanding, in 
its rest system, at a speed close to that of sound 
in a relativistic plasma (v s =c/yo). CBs contin- 
uously intercept electrons and nuclei of the ISM, 
ionized by the GRB's 7-rays. Rapidly, such an 
expanding CB becomes 'collisionless': its radius 
becomes smaller than the interaction length be- 
tween the constituents of the CB and the ISM. 
But a CB still interacts with the charged ISM 
particles, for it contains a strong magnetic field. 

Assume that the ISM particles entering a CB's 
magnetic mesh are trapped, slowly re-exiting by 
diffusion. Then, the CB's mass increases as: 

M CB « Af 7o/(/3 7 ), p = vV - 1/7, (2) 

and, for an approximately hydrogenic ISM of lo- 
cal density drii n , the CB's LF decreases as: 

d 7 /(/3 3 7 3 ) « - (m p /M 7o) dn in ( 7 ). (3) 

We approximate a CB, in its rest system, by a 
sphere of radius #(7). The ISM particles that are 
intercepted — isotropized in the CB's inner mag- 
netic mesh, and re-emitted — exert an inwards 
force on it that counteracts its expansion. The 
computed behaviour of R(j) is well described by: 

^cb(7) « [7o/(/?7)] 2/3 , Ro ~ 10 14 cm. (4) 

The distances before a CB stops and blows up 
range from a fraction of a kpc to many kpc, de- 
pending on the density profile they encounter. 

The interactions of a CB and the ISM consti- 
tute the merger of two plasmas at a large relative 



LF. This merger is very efficient in creating tur- 
bulent currents and magnetic fields (MFs) within 
the CB. We assume that the MFs, as the CB 
reaches a quasi-stable radius, are in 'equiparti- 
tion': the MF energy density equals that of the 
ISM particles the CB has temporarily phagocy- 
tized. This results in a MF |TT]: 

B CB = 3 G (7/IO 3 ) [n in /(10- 3 cm" 3 )] 1/2 , (5) 

where n m is normalized to a typical value in the 
'superbubbles' in which most SNe and GRBs are 
born. This result for .8(7) is supported by the 
analysis of the spectral evolution of GRB AGs. 

Charged particles interacting with turbulcntly 
moving MFs tend to gain energy: a 'Fermi' ac- 
celeration process. A 'first-principle' numeri- 
cal analysis [22 of the merging of two plasmas 
demonstrates the generation of such chaotic MFs, 
and the acceleration of a small fraction of the in- 
jected particles to the approximate spectrum: 

dN/dj A cx 7 ; 2 - 2 ©(7^ -7)©[7max(7)-7j. 

7ma x(7) - 10 5 7o /3 i z / A ) 7 1/3 , (6) 

with Z and A the nuclear charge and mass. The 
first 9 function reflects the fact that it is much 
more likely for the light particles to gain than 
to lose energy in their elastic collisions with the 
heavy 'particles' (the CB's turbulent MF do- 
mains). The second is the Larmor cutoff im- 
plied by the finite radius and MF of a CB. But for 
the small dependence on the nuclear identity (the 
factor Z/A), the spectrum of Eq. JBJ is universal. 

The average number of pulses in a GRB's 7- 
rays is ~ 6. Thus, the total energy of the two jets 
of CBs emitted by a core-collapse SN is: 

Eftets] ~ 12 70 N B m p c 2 ~ 2xl0 51 erg. (7) 

Practically all of this energy will, in our model, 
be transferred to CRs and the MFs they produce. 

Let i?sN ~ 2 per century be the SN rate in our 
Galaxy. In a steady state, if the low-energy rays 
dominating the CR luminosity are chiefly Galac- 
tic in origin, their accelerators must compensate 
for the escape of CRs from the Galaxy. The Milky 
Way's luminosity in CRs is then: 

LcRW-RsN^|jets]wl.3x lO^ergs" 1 . (8) 
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3. Relativistic magnetic rackets 

The essence of our theory of CRs is kinemati- 
cal and trivial. In an elastic interaction of a CB 
at rest with ISM electrons or ions of LF 7, the 
light recoiling particles (of mass m) retain their 
incoming energy. Viewed in the ISM rest system, 
they have a spectrum extending, for large 7, up 
to E ~ 2j 2 mc 2 . A moving CB is a gorgeous 
Lorentz-boost accelerator: the particles it elasti- 
cally scatters with ~ 100% efficiency reach up to, 
for 7 = 70 = (1 to 1.5)xl0 3 , an A-dependent knee 
energy E kncc (A) ~ (2 to 4) x 10 15 AeV. 

A particle of LF 7 entering a CB at rest may be 
accelerated by elastic interactions with the CB's 
turbulent plasma. Viewed in the rest system of 
the bulk of the CB, the interaction is inelastic and 
the particle may re-exit with a LF up to 7max ~ 
10 7 7; see Eq. (@J. Boosted by the CB's motion 
the spectrum of the scattered particles extends to 
7 CR ~2x 10 7 7 2 , in the ultra-high-energy (UHE) 
domain, for 7^ 70 ~ 10 3 . This powerful Fermi- 
Lorentz accelerator completes our theory of CRs. 

The calculation of the CR spectra takes 2.6 
pages 0, which I do not have here. For 7>2, to 
a good approximation, the elastic contribution to 
the CR flux of a nucleus of ISM abundance n. is: 
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where /3 C onf is the same 'confinement' index as in 
Eq. H10[) : dFeiast / dj cn depends on the priors n A , 
/Sconf > and 70, but not on any parameter specific to 
the mechanism of CR acceleration. The inelastic 
contribution is equally simple [7|. 

The source spectrum of a CR nucleus is the sum 
of an elastic and an inelastic flux, illustrated, for 
protons, in Fig. EI The elastic flux is larger than 
the inelastic one by a factor / ~ 10 at the nomi- 
nal position of the proton's knee. This ratio / is 
the only input for which we have no 'prior' infor- 
mation, and the only parameter to choose in an 
unprcdetermined range. The other parameter in 
Fig. |21 N p) is the norm of the proton inelastic flux 
at the proton's knee. Albeit within large errors, 
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Figure 2. Contributions to the E 2 -weighed pro- 
ton source spectrum. Notice the predicted 'hip'. 

N p (or N p /) is determined from the luminosity 
of Eq. (JSJ, and the Galactic (or universal) rate 
of SNe. In the domains wherein they behave as 
power laws, the elastic (inelastic) source spectra 
have indices /3 e i ast = 13/6 (~/3 clast + 0.3). 

4. Tribulations of a Cosmic Ray 

The source and local fluxes of the main CR 
elements are shown in Fig. |21 These fluxes differ 
because, on its way from its source to the Earth, a 
CR is influenced by the ambient magnetic fields, 
radiation and matter. An extragalactic CR is also 
affected by cosmological redshift. Three types of 
CR 'tribulations' must be considered: 

Interactions with magnetic fields. Fluxes of 
CRs of Galactic origin, below the free-escape (an- 
kle) energy, E ank i e (Z) ~ Z x (3xl0 18 eV), are en- 
hanced proportionally to their confinement time: 

r conf cx [ZGeV/{cp)f— f , /3 conf ~ 0.6 ± 0.1. (10) 

At higher energies CRs escape or enter the Galaxy 
practically unhindered |23| . Lower-energy extra- 
galactic CRs entering the Galaxy must overcome 
the effect of its exuding magnetic wind [Jj. 

Interactions with radiation, significant for CRs 
of extragalactic origin. The best studied one 
is the 'GZK' 7r-photoproduction on the CMB. 
Pair (e + e~) production is analogous. Photo- 
dissociation on the infrared CBR is also relevant. 

Interactions with the ISM are well understood. 
CR spallation gives rise to 'secondary' CRs. 
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Figure 3. Predicted spectra for the abundant 
elements and groups. The vertical scales are 
E 3 dF/dE. (a): The source spectra, with a com- 
mon arbitrary normalization, (b): The CR spec- 
tra at the location of the Earth. Notice that both 
the horizontal and vertical scales arc different. 

5. Results 

Relative abundances. It is customary to dis- 
cuss the composition of CRs at a fixed energy 
E A — 1 TeV. This energy is relativistic, below 
the corresponding knees for all A, and in the do- 
main wherein the fluxes are dominantly elastic 
and very well approximated by a power law of in- 
dex /3 t h = /?ciast+/3conf — 2.77. Expressed in terms 
of energy (E A cx Ay), and modified by the con- 
finement factor of Eq. Eq. 10 becomes: 

dF ohs /dE A cx fi A A^- 1 E A ^\ (11) 

with n A an average ISM abundance. At fixed 
energy the prediction for the CR abundances rel- 
ative to hydrogen is: X CIi (A) = (n A /n p ) A 1 - 77 . 
The results, for input n A 's in the 'superbubbles' 
wherein most SNe occur, are shown in Fig. 
Eq. (|TT)l snugly reproduces the large enhance- 
ments in the heavy CR abundances relative to 
hydrogen, with respect to solar or superbubble 
abundances (e.g. A 177 = 1242 for Fe). 




H He C N O Ne Mg Si S Ar Ca Fe Ni 
ELEMENT 

Figure 4. The relative abundances of primary CR 
nuclei, from H to Ni around 1 TeV. The stars 
are solar-ISM abundances. The circles are the 
predictions, with input superbubble abundances, 
a bit more 'metallic' than the solar ones. The 
squares are the CR observations. 




io" io 12 10" io" 10" io' 5 10" 10" io ,s io ! ° 10 



Primary Energy E (eV) 

Figure 5. The all-particle CR spectrum. The ver- 
tical scale is E 3 dF/dE. Some of the high-energy 
data disagree with others, and with our theory, 
which opts for the results of fluorescence detec- 
tors. The two curves correspond to two very dif- 
ferent choices for the 'penetrability' of the Galaxy 
to extragalactic CRs of energy below the ankle 
and have slightly different values of two priors 
(both within errors): the width of the 70 distri- 
bution and /3 C onf, defined in Eq. Q1U[I. 
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Figure 6. The i? 3 -weighed UHECR spectrum 
The colour coding is as in Fig. [21 



The all-particle spectrum is shown in Fig. [3] 
The curves in this figure (and in later ones) corre- 
spond to two very different choices of the Galaxy's 
'penetrability' to extragalactic CRs [7]. 

The UHECR all-particle spectrum is shown in 
Fig. |SJ At E = -Eankic the extragalactic contri- 
bution is ~ 1/2 of the data. The flux shape and 
norm, at this energy or above it, are approximate 
but 'absolute': they are the 'look-back-time' inte- 
gral of the CR flux due to SNe in other galaxies. 
The shape of the high-energy end-point and the 
height of the hump reflect not only the GZK cut- 
off, but also the 'Larmor' cutoff for protons. 

The knee region. There are recent KASKADE 
data attempting to disentangle the spectra of in- 
dividual elements or groups in the knee region. 
Our predictions for the spectra of H, He and Fe 
are shown in Fig. [7| The green line in the pro- 
ton entry has a narrow 70 distribution, the red 
and blue lines have wider distributions (all within 
the errors in this 'prior') and correspond to the 
equally coloured lines in Figs. El and EJ 

The low-energy spectra. In Fig. [Sjwe show the 
weighted spectra E^ 5 dF/dEk of protons and a 
particles, as functions of Ek, the kinetic energy 
per nucleon. The data were taken at various times 
in the 11-year solar cycle. The most intense fluxes 
correspond to data near a solar minimum. The 
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Figure 7. Measurements of individual-element 
CR spectra in the 'knee' region [251 - T ne ver ~ 
tical scales are E 2 5 dF/dE. Top: protons; mid- 
dle: a particles; bottom: iron nuclei. The data 
were kindly provided to us by K.H. Kampert. The 
colour coding is as in Figs.|3]and[n] The green line 
in the proton entry corresponds to a narrow prior 
70-distribution. The a data and, more so, the Fe 
data, show the predicted low-energy 'hips', see 
also Figs. [5] and El Notice the dependence of the 
KASKADE results on the chosen 'shower Mon- 
teCarlo'. The data are insufficient to tell apart 
-Eknec oc A (our prediction) from £"kneo oc Z (some- 
times quoted as the standard-theory result, but 
correct iff the limit is a Larmor cutoff and CR 
nuclei were accelerated up to their knees). 
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Figure 8. The very-low-energy fluxes of protons 
and a particles at various times in a solar cycle. 
The f 998 data are close to solar-minimum time. 
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Figure 9. Mean logarithmic mass of CRs. Data 
points were compiled by Hoerandel 25 from 
experiments measuring electrons, muons, and 
hadrons at ground level. The colour-coded lines 
correspond to the same choices of priors as in 
Figs. EH3 The compilation of data was kindly 
provided to us by K.H. Kampert. 



curves do not model the effects of the solar wind. 
They should agree best with the solar-minimum 
data, as they do, particularly for protons. The 
theoretical source spectra, dominated by the elas- 
tic contribution, are given by Eq. @ . These data 
are well below the elastic cutoff at 7 CR ~ 2 7q , and 
the predictions are independent of the 70 distri- 
bution. Thus, the shape of the theoretical spectra 
is, in this energy domain, parameter- free. 

Rough measures of CR composition. The evo- 
lution of the CR composition with energy is of- 
ten presented in terms of the mean logarithmic 
atomic weight (In A) , or of the depth into the at- 
mosphere of the 'maximum' of the CR-generated 
shower, X max . The predicted (In A(E)} is com- 
pared with relatively low-energy data in Fig. El 
Results for X max (E), arc shown in Fig.HT)! 

The predicted (In A(E)) at all energies, shown 
in Fig. 1111 shows how, at very high energies, 
the flux is once more Fe-dominated: lighter el- 
ements have reached their GZK, acceleration and 
Galactic-escape cutoffs. Naturally, this predic- 
tion is sensitive to the details of Galactic escape 
and extragalactic photodissociation. 




Figure 10. The depth of shower maximum as a 
function of energy. The data are from a compi- 
lation in Ref. The colour-coded lines cor- 
respond to the same choices as in Figs. EH3 an d 
OH The theoretical lines are constructed with the 
simplified method of Wigmans 27. , arguably as 
good as any currently-used MonteCarlo. 
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Figure 11. Predicted (\nA(E)) at all energies. 



Cosmic-ray electrons. Electrons and nuclei are 
accelerated by the 'magnetic-racket' CBs in the 
same manner |7I13| . The functional form of their 
source spectra is therefore the same, dF s /dj cx 
7 /3ola3t , in the range 10 < 7 < 10 6 . Electrons 
lose energy much more efficiently than protons 
in their interactions with radiation, MFs and the 
ISM. The rates -dE/dt cx E a of their various 
mechanisms of energy loss have different a's. For 
Coulomb losses a — 0; for bremsstrahlung a — 1; 
for inverse Compton scattering (ICS) and syn- 
chrotron losses at the relevant energies, a = 2. 
At sufficiently high energy, the radiative energy 
loss dominates the others. In this domain, the 
steady-state solution of the equation describing 
the radiation-modified electron spectrum, for a 
source dF s /dE cx E-P'*"*, is simply |28I13| : 

dFjdEcxE-l 3 *; e = /3 elast + 1 « 3.17. (12) 

This result agrees with the observed slope of the 
CRE spectrum; see Fig.^1 The best- fitted value 
above GeV is /Us = 3.2 ± 0.10, and the fit 
is excellent if the experiments are recalibrated to 
the same flux at high energy |13| . 

The GBR. The existence of a diffuse gamma 
background radiation, suggested by data from the 
SAS 2 satellite, was confirmed by the EGRET 
instrument on the Compton Gamma Ray Ob- 
servatory We call 'the GBR' the diffuse 
emission observed by EGRET by masking the 
galactic plane at latitudes |b| < 10°, as well as 
the galactic centre at |b| < 30° for longitudes 
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Figure 12. The CRE spectrum, compiled in |2S|. 
The line is the central result of a power fit to the 
higher-energy data; its slope is 3.2±0.1 I3TJI • The 
different slope below ~ 6 GeV results from the 
diffusion of electrons in the Galactic MF 0. 



|1| < 40°, and by extrapolating to zero column 
density, to eliminate the 7r° and bremsstrahlung 
contributions to the observations and to tame the 
model-dependence of the results. Outside this 
'mask', the GBR flux integrated over all direc- 
tions, shown in Fig. 1131 is well described by a 
power law dF 7 /dEocE~ 2 - w±a - 03 [30]. 

The EGRET GBR data show a significant devi- 
ation from isotropy, correlated with the structure 
of the Galaxy and our position relative to its cen- 
tre [2H1- Contrariwise, the GBR's spectral index 
is uncannily directionally uniform. These facts 
suggest a GBR that is partially local, as opposed 
to dominantly cosmological, and a common origin 
for the Galactic and extragalactic contributions. 

In |28ll3j we have analyzed the directional and 
spectral properties of the EGRET data and con- 
cluded that the GBR is produced by ICS of CREs 
on starlight and the CBR. It has comparable con- 
tributions from CREs in a Galactic halo of di- 
mensions akin to the hadronic-CR confinement 
volume (a directional and local source) and from 
other galaxies (an isotropic cosmological compo- 
nent). Thus, the GBR is a CR 'secondary'. Its 
spectral index is the same for the local and cos- 
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Figure 13. The GBR spectrum, measured by 
EGRET 30 . The line is the central result of 
a power-law fit of slope 2.10 ± 0.03. 

mological contributions. 

If produced by ICS by electrons with the spec- 
trum of Eq. H12(l . the GBR has a spectrum: 

dF*/dEocE-^, /3 7 = (/3 e -l)/2~2.08. (13) 

The predicted photon spectral index |28I13| coin- 
cides with the measured one, 2.10 ± 0.03 [3Tj|. 

Other predictions. Our CR theory explains 
other observations that I have no space to discuss 
here: the slight differences between the slopes of 
the spectra of the hadronic CRs, the deduced con- 
finement time and volume of CR electrons and 
nuclei in the Galaxy [7j, and the normalization 
and directional dependence of the GBR flux ^21 ■ 

6. Discussion and conclusions 

Our theory is incomplete in several respects. 
The ejection of CBs in episodes of accretion onto 
a compact object is supported by observations, 
but not fundamentally understood (this is also 
the case for the ejection of SN shells). Because 
CBs deposit CRs along their kpc-long trajecto- 
ries, CR diffusion does not play a crucial role, 
but we have not studied it in detail (in the stan- 
dard theory diffusion results in directional asym- 
metries that are not observed, and calls for ad-hoc 



remedies, such as 'CR reacceleration'). The tem- 
porary confinement of CRs — in the Galaxy or 
within a CB — is not fully understood. Neither is 
the 'penetrability' of the Galaxy to extragalactic 
CRs of energy below the ankle. We have studied 
the dynamics of the expansion of CBs, but not 
modelled it in minute detail. We contend that a 
good fraction of the original energy of CRs ends 
up in the production of 'equipartition' MFs |31|. 
but we cannot predict this, nor determine the ef- 
fect that it may have on the CR spectral shape. 

In spite of the above limitations, we have 
demonstrated how our simple and single accelera- 
tors — cannonballs — are effective at all observed 
energies. The mechanisms of CR acceleration, 
particularly the 'elastic' one, could hardly be sim- 
pler: a fast massive object slows down by kicking 
out of its way the light particles it encounters. 

Most of our results are 'robust' in that — within 
very large brackets — they do not depend on the 
specific choices of parameters and priors: 

• An all-particle piecewise power-law spectrum 
with four features: two steepenings at the knee 
and the second knee, a softening at the ankle, 
and an end-point at the roughly-coincident GZK 
and proton-acceleration cutoffs. 

• An UHECR flux above the ankle, which is 
predicted — to within a factor of a few — and 
otherwise parameter-free. 

• A composition dependence at 1 TeV with the 
observed trend, so different from that of the ISM. 

• A very low-energy flux whose spectral shape is 
independent of any CB-model 'prior' parameters. 

• Individual-element knees that scale like A and 
occur at the predicted energies. 

• A non-trivial shape of the individual knees: an 
abrupt decrease in flux, followed by a spectrum 
steeper than that below the knee. 

• An ankle with the observed shape. The dom- 
inantly Galactic-Fe flux below it and the dom- 
inantly extragalactic-proton flux above it are 
comparable in magnitude at the estimated es- 
cape 'ankle' energy [231 of Galactic protons. 

• A composition dependence that is almost 
energy-independent below the knee becomes 
'heavier' from the knee to the second knee, 
'lighter' again above it, and finally heavier at 
yet unmeasured ultra-high energies. 
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• An 'extended' distribution of CR sources along 
CB trajectories emerging from the central realms 
of the Galaxy, where most SN explosions oc- 
cur, implying a CR flux at the Earth's location 
with a much smaller and less energy-dependent 
anisotropy than that of standard models of CRs. 

• Predictions for the values of a related set of ob- 
servables: the CR luminosity, confinement time 
and volume of the Galaxy, the spectral indices of 
CR electrons and of the diffuse GBR. 

Our results describe the observed properties of 
hadronic non-solar CRs very well from the lowest 
energies to ~ 10 10 GeV. Above that energy and 
up to the highest observed energies, ~ 10 11 GeV, 
our theory opts for the data gathered with fluo- 
rescence detectors. Overall, the energy range for 
which the theory is successful covers ten decades 
and the flux extends over three times as many. 
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